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In this paper the single and multibunch omittance growths in the main linac of a linear 
collider are analytically treated in analogy to the Brownian motion of a molecule, and the 
. analytical formulae for the emittance growth due to accelerating structure misalignment 

errors are obtained by solving Langevin equation. The same set of formulae is derived also 
by solving directly Fokker-Planck equation. Comparisons with numerical simulation results 
' are made and the agreement is quite well. 
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To achieve the required luminosity in a future e + e linear collider one has to produce two 
colliding beams at the interaction point (IP) with extremely small transverse beam dimensions. 
According to the linear collider design principles described in ref. 1, the normalized beam 
emittance in the vertical plane (the normalized beam emittance in the horizontal plane is larger) 
1— i ' at IP can be expressed as: 

_ n*r e 

1€y ~ 374<%a 4 1 j 

where 7 is the normalized beam energy, r e = 2.82 x 10 -15 m is the classical electron radius, 
a = 1/137 is the fine structure constant, 5g is the maximum tolerable beamstrahlung en- 
ergy spread, and n 7 is the mean number of beamstrahlung photons per electron at IP. Taking 
5* B = 0.03 and n 7 = 1, one finds "fe y = 8.86 x 10 _8 mrad. To produce beams of this small trans- 
verse emittance damping rings seem to be the only known facilities which have the potential to 
do this job. The questions now are that once a beam of this small emittance is produced at the 
exit of the damping ring, how about the emittance growth during the long passage through the 
accelerating structures and the focusing channels from the damping ring at the beam energy of 
few GeV to the IP with the beam energy of a few hundred of GeV, and how to preserve it ? 
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Many works have been dedicated to answer these questions [|[ H [f| f| • To start with, in 
sections 2, 3, and 4, we consider the short range wakefield induced single bunch emittance 
growth and try to calculate the emittance growth by using two different methods, and show 
that the two methods give the same results. Since the number of accelerating structures in the 
main linac of a linear collider is very large, the transverse random kicks on the particles can be 
described statistically. Firstly, we make use of the analogy between the transverse motion of 
a particle in a linear accelerator with the Brownian motion of a molecule, which are governed 
by Langevin equation. Secondly, we solve directly Fokker-Planck equation. What should be 
noted is that both methods are physically consistent. As a natural extension, multibunch case 
is treated in section 5. Comparisons with some numerical simulation results are made in section 
6. 



2 Equation of transverse motion 

The differential equation of the transverse motion of a bunch with zero transverse dimension is 
given as: 

d 2 y(s,z) 1 dj(s,z) dy(s,z) 2 

, 2 + r — a a + fc ( s ' z > y ( s ' z > 

ds L 7{s,z) as as 

1 f°° 

e 2 N e p(z')W±(s,z' -z)y(s,z')dz' (2) 



m c 2 7(s 

where k(s, z) is the instantaneous betatron wave number at position s, z denotes the particle 
longitudinal position inside the bunch, and p(z')dz' = 1. Now we rewrite eq. [2] as follows: 

d 2 y(s,z) dy{s,z) 2 

— ja — + r — j + z) y(s, z) = A (3) 

as z ds 

where V = 7 1°^ , G = — e ^ z , n \ , E z is the effective accelerating gradient in the linac, 

A = e2Ne Zc^(p) S,0) > W ± (s,z) = f™p(z')W ± (s,z' - z)dz', and y(s,0) is the deviation of the 
bunch head with respect to accelerating structures center. In this section we consider the case 
where the injected bunch, quadrupoles and beam position monitors are well aligned, while the 
accelerating structures are misaligned. As a consequence, y(s, 0) is a random variable exactly the 
same as random accelerating structure misalignment errors with < y(s,0) >= (< > denotes 
the average over s). If we take z as a parameter and regard T, k(s,z), and A as adiabatical 
variables with respect to s, eq. ||] can be regarded as Langevin equation which governs the 
Brownian motion of a molecule. 



3 Method one: Langevin equation 

To make an analogy between the movement of the transverse motion of an electron and that of 
a molecule, we define P = e ^^^^^y 2 > an d regard y(s, 0)P as the particle's "velocity" random 

increment (A^) over the distance l s , where l s is the accelerating structure length. What we 
are interested is to assume that the random accelerating structure misalignment error follows 
Gaussian distribution: 

^•°^vk^(- v W) (4) 

and the velocity (u) distribution of the molecule follows Maxwellian distribution: 

I ffl / 711 'U^ \ 
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where m is the molecule's mass, k is the Boltzmann constant, and T is the absolute temperature. 
The fact that the molecule's velocity follows Maxwellian distribution permits us to get the 
distribution function for Al s ||: 

where 

kT 

q = T— 7 

m 



By comparing eq. ^ with eq. ^j, one gets: 

p2 



2a v = -nST ( 8 ) 



or 

m 2l s T K ' 

Till now one can use all the analytical solutions concerning the random motion of a molecule 
governed by eq. || by a simple substitution described in eq. ||. Under the condition, k 2 (s, z) >> 
^- (adiabatic condition), one gets ||: 

< 2/ 2 >= — TJT — \ + f^o TJT — \1 fcos(A;is) + sm(kis)) exp(-rs) 



+ 



mk 2 (s,z) V mk 2 (s,z)J \ 2k\ 
cr 2 y l s f e 2 N e W ± (z)\ 2 



27(5, z)^f(0)Gk 2 (s, z) \ rriQC 
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Vl ~ „ , s4^TJ7 ^ ( e2jV ^ (Z) l I ( «»(M + ^ S m(fci S )fexp(-r S ) (10) 



27(s,z)7(0)G/c 2 ( s ,z) ^ m c 2 / / V 2fc 
kT Jxii^s ^ z^j f 2 kT 



ggjg / e 2 iW±(^) N 2 



+ 



2 7 (s,z)7(0)GA; 2 (s,z) ^ m c 2 

^ (* - ^JTO I I -"<M«P<-r.) (") 



2^ 



/ 



^ fcl V^I^S, z) 2 



y 2 ,^ fcos(A;is) + ^-sin(/cis)^ exp(-rs) 



k{s,z) 2 ( a 2 l s f e 2 N e W ± (z) 

fa I 2^(s,z)j(0)Gk 2 (s,z) I m c 2 



2 



cos(fcis) H — — sin(fcis) ) exp(— rs) (12) 



2 A; 



where fa = ^k{s 1 z) 2 — \T 2 . The asymptotical values for < y 2 >, < y' 2 >, and < yy' > as 
s — > 00 are approximately expressed as: 

,,2._ ^ _ gg, / e 2 iV e ^ ( z) \ 2 

mk 2 (s,z) 2j(s, z)j(0)Gk 2 (s, z) \ tuqc 2 J 



3 



< y' 2 >= k A (s,z) < y l >- 



e 2 N e W ± (z) 



2 7 (s,z)7(0)G I m c 2 



< yy' >= 



(14) 
(15) 



Inserting eqs. 13, 14, and 15 into the definitions of the r.m.s. emittance and the normalized 



1/2 



r.m.s. emittance shown in eqs. 16 and 17: 

e T ms = (< y 2 >< y' 2 > - < yy' > 

e n ,rms = 70, z) (< y 2 >< y' 2 > - < yy' > 



A 1/2 



one gets 



and 



e 2 N e W ± (z) 



2j(s, z) , y(0)Gk(s, z) \ vtiqc 2 



e 2 N e W ± (z) 



2j(0)Gk(s,z) I m c 2 



(16) 
(17) 

(18) 
(19) 



The effects of energy dispersion within the bunch can be discussed through 7(5, z), k 2 (s, z), 
and W±(z), such as BNS damping ||. From eqs. 13, 18, and 1£ one finds that there are 
three convenient types of variations of k(s, z) with respect to s. If one takes k 2 (s, z) r f(s, z) = 
k 2 (0, 2)7(0, z), one gets < y 2 > independent of s. If one takes, however, k(s, z)^f(s, z) = 
k(0, 2)7(0, z), e r ms is independent of s, and finally, if k(s,z) = k(Q,z), one has tn,rms is in- 
dependent of s. One takes usually the first scaling law in accordance with BNS damping. To 
calculate the emittance growth of the whole bunch one has to make an appropriate average over 
the bunch, say Gaussian as assumed above, as follows: 



bunch = Ho P(z')en,rms(z')dz' 

IZoP(z')dz' 



(20) 



To make a rough estimation one can replace p{z) by a delta function 6(z — z c ), and in this case 
the bunch emittance can be still expressed by eq. [l9| with W±(z) replaced by W±(z c ), where z c 
is the center of the bunch. 



4 Method two: Fokker-Planck equation 

Keeping the physical picture described above in mind, one can start directly with Fokker-Planck 
equation which governs the distribution function of the Markov random variable, y'\ 

S -EM1 = -^ ( AF(^)) + \^ { DF M) ) (21) 

with 

A = K< y >> (22) 

D= «^f» m 

where Ay' is the increment of y' over l s , and << >> denotes the average over a large number of a 
given type of possible structure misalignment error distributions (in numerical simulations, this 
average corresponds to the average over the results obtained from a large number of different 
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seeds, for a given type of structure misalignment error distribution function, say, Gaussian 
distribution). From eq. one gets the increment of y' over l s : 



Ay' = (l-exp(-^)y + n s 



In consequence, one obtains: 



«(A 2/ ')»~(l-exp(-^))y / 

« (Ay') 2 »« (1 - exp(-^))V 2 + « (M s ) 2 » exp(-I7 s ) 
where << Tl s >>= has been used. Inserting eqs. ^ and|3^ into eq. one gets: 



I 



dF(s,y>) n a/fM 

^— = -(l-exp(- T )) - 



>l-eM- I t)) 2 d 2 (y' 2 F(s,y')) , « (AZ S ) 2 » 



+ 



2 dy' 2 2 

Multiplying both sides with y' 2 and integrating over y', one has: 

d < y' 2 > 



exp(-rZ s ) 



d 2 F 

dy' 2 



(1 - exp(-rZ s )) < y' 2 > + « (Al s )' 2 » exp(-r/ s 



ds 



Assuming that Tl s << 1, eq. 28 is reduced to: 



h d< f > = -rh <y' 2 > + « (Ah) 2 » 

ds 



Solving eq. 29, one gets: 



>2 



< y > 



« A 2 » l s 
2T 



(l-exp(-rs)) + exp(-ra)^ 
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where y' is the initial condition. Apparently, when s — > oo, one has: 



(24) 

(25) 
(26) 



(27) 



(28) 



(29) 



(30) 



< y' 2 >oo= 



<< A 2 >> l s 
2T 



2 yl S 



e 2 N e W ± (z) 



2 7 (s,z)7(0)G V m c 2 



(31) 



where a 2 =<< y(s,0) 2 ». Eq. 31 is the same as what we have obtained in eq. 14. In fact, by 
solving directly Fokker-Planck equation we obtain the same set of asymptotical formulae derived 
in section 3. 



5 Multibunch emittance growth 



Physically, the multibunch emittance growth is quite similar to that of the single bunch case, 
and each assumed point like bunch in the train can be regarded as a slice in the previously 
described single bunch. Obviously, the slice emittance expressed in eq. 19 should be a good 



starting point for us to estimate the emittance growth of the whole bunch train. Before making 



use of eq. |19| let's first look at the differential equation which governs the transverse motions of 
the bunch train: 



d_ 

ds 



ln(s) 



dy n 
ds 



+ 7r 



(s)k 2 n y r , 



3 2 iV f 



n-1 



W T ((n - i)s b ) yi 



(32) 
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where the subscript n denotes the bunch number, s& is the distance between two adjacent 
bunches, N e is the particle number in each bunch, Wt(s) is the long range wakefield produced 
by each point like bunch at distance of s. Clearly, the behaviour of the ith bunch suffers from 
influences coming from all the bunches before it, and we will treat one bunch after another in an 
iterative way. First of all, we discuss about the long range wakefields. Due to the decoherence 
effect in the long range wakefield only the first dipole mode will be considered. For a constant 
impedance structure as shown in Fig. 1, one has: 

where a z is the rms bunch length (a z is used to calculate the transverse wake potential, and 
the point charge assumption is still valid), u% and Q\ are the angular frequency and the loaded 
quality factor of the dipole mode, respectively. The loss factor k\ in eq. |33] can be calculated 
analytically as 0: 

hJ 2 (^-a) , x9 
e TvDR 2 J$(uii) 

S(x) = ™£ (35) 
x 

xi = -m (36) 

where R is the cavity radius, a is the iris radius, h is the cavity hight as shown 
in Fig. [l], and un = 3.832 is the first root of the first order Bessel function. 




Figure 1: A disk-loaded accelerating structure. 

To reduce the long range wakefield one can detune and damp the concerned dipole mode. The 
resultant long range wakefield of the detuned and damped structure (DDS) can be expressed as: 



W t ,dds(s) 



1 ^ 2ck 1:i s ( uu 



Nr. ^ Wl,iO? 



exp 



2c 2 



(37) 



where N c is the number of the cavities in the structure. When iV c is very large one can use 
following formulae to describe ideal uniform and Gaussian detuning structures ||: 
1) Uniform detuning: 



W T ,i,u = 2<K>sm[ — — exp - 

V c J (vrsA/i/c) V < Q >i c 



(38) 



2) Gaussian detuning: 



W TXG = 2 < K > sin ( 2 *<f >s ) e -2.^, - < . 



exp 



7T < / > S 

<Q>ic 



(39) 
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where K = ■pfa, /i = g£, A/i is full range the synchronous frequency spread due to the 
detuning effect, aj is the rms width in Gaussian frequency distribution. 

Once the long range wakefield is known one can use eq. 13 to estimate < yf > in an iterative 
way, and the emittance of the whole bunch can be calculated accordingly as we will show later. 
For example, if a bunch train is injected on axis {y n = 0) into the main linac of a linear collider 
with structure rms misalignment a y , at the end of the linac one has: 

,2 



< yf >= o (40) 



2 7 ( S )7(0)GA:2( s )( moC 2)2 > 



{al + ^)e 2 N e \W T (2s b )\ + j (a* + S § 2l )e 2 N e \W T (s b ) 



2 

I, 



<%2> ' 2 7 ( S ) 7 (0)GA;2( s )(moc2)2 (J2) 



and in a general way, one has: 



2 

's 



Epj + I < ~W>)e 2 N e \W T {(i - j)s h )\) I 
2 1 {s) 1 {V)Gkl{s)(m 6 



Finally, one can use the following formula to estimate the projected emittance of the bunch 
train: 



train _ l(. s )^( s , 



e;;™ = • ' Ar • ' J2 < y 2 i > ( 44 ) 



where k(s) = k n (s) (since the bunch to bunch energy spread has been ignored), and k(s) is the 
average over the linac. 

It is high time now for us to point out that the analytical expressions for the single and 
multibunch emittance growths established above give the statistical results of infinite number of 
machines with Gaussian structure misalignment error distribution, which corresponds to using 
infinite seeds in numerical simulations. 



6 Comparison with numerical simulation results 

To start with, we take the single bunch emittance growth in the main linac of SBLC §1 for 
example. The short range wakefields in the accelerating S-band structures are obtained by 
using the analytical formulae Q and shown in Fig. |2[ In the main linac the beam is injected at 
3 GeV and accelerated to 250 GeV with an accelerating gradient of 17 MV/m. The accelerating 
structure length l s =6 m, the average beta function (3(s) is about 70 m (k(s,z) = for 
smooth focusing), the bunch population N e = 1.1 x 10 10 , the bunch length a z = 300 fim, and 
and corresponding dipole mode short range wakefield W±(z c ) = 338 V/pC/m 2 . Inserting these 
parameters into eq. [H], one finds e n>rms = 8.66 x a 2 z . If accelerating structure misalignment 
error a y = 100 ^m, one gets a normalized emittance growth of 8.66xl0 -8 mrad, i.e., 35% 
increase compared with the designed normalized emittance of 2.5 x 10 -7 mrad. The analytical 
result agrees quite well with that obtained from numerical simulations |J. Now, we apply the 
analytical formulae established for the multibunch emittance growth to SBLC, TESLA and NLC 
linear collider projects where enormous numerical simulations have been done. The machine 
parameters are given in Tables 1 to 4 which have been used in the analytical calculation in this 
paper. Firstly, we look at SBLC. 
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Figure 2: The short range wakefields of SBLC type structure with a z = 300 /iin, and the 
beam iris a = 0.0133 m. (a) and (b) are the bunch current distributions, (c) is monopole the 
longitudinal wakefield. (d) is the dipole transverse wakefield at r = a. 

Fig. H shows the "kick factor" K defined in eqs. 38 and 39 vs the dipole mode frequency. 



Figure 3: The K{ vs dipole mode frequency (SBLC). 



Fig. |4|(a) gives the long range transverse wakefield produced by the first bunch at the locations 
where find the succeeding bunches, while Fig. ||(b) illustrates the square of the rms deviation 
of each bunch at the end of the linac with the dipole loaded quality factor Q\ = 2000. The 
corresponding results for Q\ = 10000 are shown in Fig. [B]. The normalized emittance growths 
compared with the design value at the interaction point {t^rms ,IP = 2-5 x 10 -7 mrad) are 32% 
and 388% corresponding to the two cases, respectively as shown in Table 4, which agree well 
with the numerical results Q. 



S 



Figure 4: (a) the long range dipole mode wakefield vs the number of bunch, (b) the yf at the 
end of linac vs the number of bunch (SBLC, Qi = 2000, a y = 100 /urn). 



Figure 5: (a) the long range dipole mode wakefield vs the number of bunch, (b) the yf at the 
end of linac vs the number of bunch (SBLC, Q\ = 10000, a y = 100 |im). 



To demonstrate the necessity of detuning cavities we show the violent bunch train blow up if 
constant impedance structures are used in spite of Q\ being loaded to 2000 as shown in Fig. ||. 
Secondly, TESLA (the version appeared in ref. 10) is investigated. From Fig. [7| one agrees that 
it is a no detuning case. 



Figure 6: (a) the long range dipole mode wakefield vs the number of bunch, (b) the yf at the 
end of linac vs the number of bunch (SBLC no detuning, Q\ = 2000, <r y = 100 /xm). 



Figure 7: The Ki vs dipole mode frequency (TESLA). 
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From the results shown in Fig. |8| and Table 4 one finds that taking structure misalignment 
error a y = 500 /xm and Q\ = 7000 one gets an normalized emittance growth of 24% which is a 
very reasonable result compared what has been found numerically in ref. 10. Thirdly, we look 
at NLC X-band main linac. To facilitate the exercise we assume the detuning is effectuated as 
shown in Fig. ^ (in reality, NLC uses Gaussian detuning). Fig. 1C shows the analytical results 



with a y = 15 /xm and Q\ = 1000. From Table 4 one finds a normalized emittance growth of 
21%. Then, we examine NLC S-band prelinac. Assuming that the detuning of the dipole mode 



is shown in Fig. 11, one gets the multibunch transverse behaviour and the normalized emittance 



growth in Fig. 12 and Table 4. Finally, in Figs. 13 and 14 we give more information about the 
emittance growth vs Q\ in NLC X-band and S-band linacs. 



Figure 8: (a) the long range dipole mode wakefield vs the number of bunch, (b) the yf at the 
end of linac vs the number of bunch (TESLA, Ql=7000, a y = 500 fim). 



Figure 9: The Ki vs dipole mode frequency (NLC X-band linac). 
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Figure 10: (a) the long range dipole mode wakefield vs the number of bunch, (b) the yf at the 
end of linac vs the number of bunch (NLC X-band linac, Qi=1000, a y = 15 /xm). 



Figure 11: The Ki vs dipole mode frequency (NLC S-band prelinac). 
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Figure 12: (a) the long range dipole mode wakefield vs the number of bunch, (b) the yf at 
end of linac vs the number of bunch (NLC S-band prelinac, Qi=10000, a y = 50 /im). 
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Figure 13: The normalized emittance growth vs Q\ with a y = 15 /jm (NLC X-band linac) 
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Figure 14: The normalized emittance growth vs Q\ with <r y = 50/um (NLC S-band linac). 



Machine 


h ( m ) 


N c 


h (GHz) 


a (m) 


D (m) 


h (m) 


R (m) 


SBLC 


6 


180 


4.2-4.55 


0.015-0.01 


0.035 


0.0292 


0.041 


TESLA 


1 


9 


1.7 


0.035 


0.115 


0.0974 


0.095 


NLC X-band 


1.8 


206 


15-16 


0.0059-0.00414 


0.00875 


0.0073 


0.011 


NLC S-band 


4 


114 


4.2-4.55 


0.015-0.01 


0.035 


0.0292 


0.041 



Table 1: The machine parameters I. 



Machine 


N e (xlO iU ) 


s b (m) 


E z (MV/m) 




N b 


Qi 


SBLC 


1.1 


1.8 


17 


300 


333 


2000,10000 


TESLA 


3.63 


212 


25 


700 


1136 


7000 


NLC X-band 


1.1 


0.84 


50 


145 


95 


1000 


NLC S-band 


1.1 


0.84 


17 


500 


95 


10000 



Table 2: The machine parameters II. 
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Machine 


7 (0) (GeV/MeV) 


7 (GeV/MeV) 


k(s) (1/m) 


a y (fim) 


SBLC 


3/0.511 


250/0.511 


1/90 j 


100 


TESLA 


3/0.511 


250/0.511 


1/90 


500 


NLC X-band 


10/0.511 


250/0.511 


1/50 


15 


NLC S-band 


20/0.511 


10/0.5111 


1/20 


50 



Table 3: The machine parameters III. 



Machine 


c train,numeri. (■ mY ,„A\ 

e n,rms (mraaj 


~train,analy. /'vnvoH^ 
t n,rms \miau.) 


(mrad) 


SBLC 


2.3xl0~ 8 ,8.8 x 10~ Y 


8.xl0~ 8 ,9.7 x 10~ y 


2.5xl0" 7 


TESLA 


~2.5xl0~ s 


5.9xl0" s 


2.5xlO~ Y 


NLC X-band 




3xl0~ 8 


1.4xl0- 7 


NLC S-band 




1.2xl0~ s 


1.4xl0 _Y 



Table 4: The normalized train emittance growth. 

7 Conclusion 

We treat the single and multibunch emittance growths in the main linac of a linear collider 
in analogy to the Brownian motion of a molecule, and obtained the analytical expressions for 
the emittance growth due to accelerating structure misalignment errors by solving Langevin 
equation. As proved in this paper, the same set of formulae can be derived also by solving 
directly Fokker-Planck equation. Analytical results have been compared with those coming 
from the numerical simulations, such as SBLC and TESLA, and the agreement is quite well. 
As interesting applications, we give the analytical results on the estimation of the multibunch 
emittance growth in NLC X-band and S-band linacs. 
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